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Abstract 

A simple analytical model is developed to analyze and explain the complex dynamics of the 
multi-peak modulation instability spectrum observed in dispersion oscillating optical fibers [M. 
Droques et al, 37, 4832-4834 Opt. Lett., (2012)]. We provide a simple expression for the local 
parametric gain which shows that each of the multiple spectral components grows thanks to a 
quasi-phase-matching mechanism due to the periodicity of the waveguide parameters, in good 
agreement with numerical simulations and experiments. This simplified model is also successfully 
used to tailor the multi-peak modulation instability spectrum shape. These theoretical predictions 
are confirmed by experiments. 
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I. INTRODUCTION 



Modulation instability (MI) is a nonlinear process in which a weak perturbation is ex- 
ponentially amplified by an intense field. MI has been investigated in many sub-fields of 
physics, especially in optics for the relative simplicity to perform experiments, in both ho- 
mogeneous and periodical media pQ. The additional degree of freedom brought by the 
periodicity has catched the attention of many research groups and has led to many theo- 
retical and experimental works, in spatial [2j [3] and temporal Kerr media such as optical 
fibers jlHT]. This physical flexibility is of particular interest since the modulation period 
can be varied from the meter range up to tens of kilometers, leading to a large window 
of investigation. In the early nineties, the rise of optical telecommunication networks has 
led to the deployment of "natural" periodic optical fiber systems due to the alternation of 
all-optical regeneration devices and/or dispersion managed lines [H [9]. In addition to the 
fundamental interest brought by these systems, it was then necessary to understand in depth 
the origin of the characteristic spurious MI sidebands [HI E] which are highly detrimental 
for telecommunications since they are generated in the GHz range [51 Ej due to their period 
in the range of kilometers. A lot of theoretical studies have therefore been initiated in this 
context 0IEHI5]. 

Very recently, dispersion management was pushed one step further with the experimental 
demonstration of MI in continuously modulated waveguides [TBJ [T7]. From these results, the 
spurious consequence of the MI sidebands due to the periodicity can be turned into benefits 
since these works show the possibility to obtain multiple parametric gain bands in the 
THz range. It should then provide another degree of freedom for designing optical systems 
requiring broad bandwidths such as in all optical signal processing systems where there is a 
growing demand. 

While a deep theoretical study of MI in periodically tapered fibers has been reported 
very recently [TS], we propose in the present work a simplified analytical treatment allowing 
to accurately describe the MI dynamics and to tailor the overall shape of its multi-peak 
spectrum. Besides providing insight into the underlying physics, our analytical treatment 
allows to derive an expression for the local linear parametric gain for the first time. Finally, 
in order to illustrate the practical interest of this simple analysis, we report experiments 
in which the multi-peak MI spectrum has been tailored to suppress one MI sideband or to 
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favor a single strong one, in good agreement with our analytical predictions. 
II. CONTEXT 

We recently reported the first striking experimental demonstration of MI in a dispersion 
oscillating fiber (DOF) [T7j. Although details can be found in [T7j, the aim of this section 
is to briefly summarize our previous results in order to facilitate the reading of the present 
paper. Figure [T] shows a measurement of the evolution of the fiber diameter along its length. 
The outer diameter follows a sine shape with a modulation amplitude of ±7 % and a period 
Z of 10 m, which results in a dispersion modulation over the fiber length z with a quasi- 
sinusoidal shape over the wavelength range of interest here, written as 



where fa = 1-2 x 1CT 27 s 2 /m is the average second order dispersion and ft^ — 1.5 x 1CT 27 s 2 /m 
is the amplitude modulation at our pump wavelength of 1072 nm. 

Full circles in Fig.[T](b) correspond to the spontaneous MI spectrum obtained by pumping 
a 120 m-long sample of this DOF (labeled DOF^l) with 2 ns pulses with a peak power P p of 
20 W at X p = 1072 nm. The solid line represents the spectrum resulting from the numerical 
integration of the generalized nonlinear Schrodinger equation (GNLSE) seeded by noise to 
accurately reproduce experimental random initial conditions (all details and parameters are 
given in [ 1 Tj ) . Both spectra are in good agreement (except for the higher experimental noise 
floor) and show the generation of multiple MI sidebands pairs [TH] spanning over more than 
10 THz. The frequency of these parametric sidebands can be roughly estimated from a 
quasi-phase-matching relation developed in the case of an infinitely long grating (6j EJ [9] 



where k is an integer, Qk is the pulsation detuning from the pump and 7 is the average 
nonlinear coefficient of the DOF. 

The dynamics of the MI process with fiber length was investigated by cutting back the 
DOF and recording output spectra. As an illustration, Fig. [IJc) shows the evolution of the 
power of the first sideband (k = 1) versus fiber length obtained from experiments (circles) 
and numerical simulations (solid line). The dynamics observed both in experiments and 
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FIG. 1: (a) Outer diameter of the DOF versus length measured during fiber drawing, (b) Ex- 
perimental (circles) and simulated (solid line) spectra obtained for a pump power of 20 W and a 
fiber length of 120 m. (c) Evolution of the power of the first sideband {k = 1) versus fiber length 
obtained from experiments (circles) and numerical simulations (solid line). Results are from |17j . 

simulations exhibits periodic regions of deamplification, which makes the side lobe power 
oscillate along the fiber around the exponential growth (expected for a perfectly phase- 
matched process). Such a dynamics is expected from quasi-phase- matched processes, but it 
differs from the one observed in second-order nonlinear crystals in which there are no regions 
of deamplification. This particular and unusual feature will now be studied in detail with 
help of a simple and intuitive analytical model, in order to provide further insight into the 
underlying physics. 



III. ANALYTICAL MODEL 



A. Parametric gain calculation 

It is well established that the MI process can be interpreted in the spectral domain as a 
four-wave mixing (FWM) process [201 I2T] . In this frame, the parametric FWM gain spec- 
trum can be obtained by studying the stability of the steady state solution against weak 
perturbations through a so-called linear stability analysis. In dispersion managed optical 
systems, this tool has allowed to analytically predict the complex multi-peak gain spec- 
trum P, [TTJ , but such an analysis does not provide any clear insight into the dynamics of 
the process nor any details about the fine evolution of the field over a single modulation 
period of the fiber. To this aim, we propose here a more intuitive explanation of the results 
from Ref. [T7] by revisiting a simplified truncated three-wave model usually aimed at de- 
scribing Fermi-Pasta-Ulam recurrence and fiber-optic parametric amplification [22T - T21] . This 
model allows to account for the relative phase variations between pump, signal and idler 
waves during propagation. In our work, it will be induced by the longitudinal variations of 
dispersion rather than pump depletion. Our starting point is the four coupled differential 
equations given by Eqs. 3 in Ref. [23]. We neglect fiber loss, we assume that the pump 
remains undepleted and that signal and idler powers, P s and Pi, are much less than the 
pump power P p over the whole DOF length. It is then easy to show that this system reduces 
to the following equations 



dFa ^ ^ = 2rfP p y/P a (Sl, z)Pi{VL, z) sin 0(fi 
dPl< f' ^ = 2 1 P pV /p i {Q, z)P s {Q, z) sin#(ft 



de(n,z) q2 



dz 



/3 2 + f3o sin 



2nz 



(3a) 
(3b) 

z j +2 1 P p (l + cos[9(n,z)]) (3c) 

where Q is the shift of the signal and idler pulsations from the pump, and 8(Q, z) describes 
the longitudinal evolution of the relative phase difference between all these waves [22] ■ The 
discrepancy between solutions of Eq. [2] and experimental/numerical values mentioned in 
Ref. [T7] can now be understood from Eq. 3c Indeed, Eq. [2] assumes that the nonlinear 



phase mismatch can be approximated by 2jP p [HI El E] , while Eq. [3c] shows that it is in fact 
equal to 27P p (l + cos [9(ft, z)] ). This does not impact the validity of the present results 
since this additional term remain low for the pump powers involved in the present study. In 



order to obtain a simple analytic solution of the set of Eqs. |3j we thus neglect the last term 



cos z)\ in Eq. 3c It physically means that we assume that the longitudinal evolution of 
the nonlinear phase mismatch term is weak as compared to the linear and uniform nonlinear 
phase mismatch terms, which is valid for low pump powers. By integrating the set of Eqs. |3j 
we find that the total accumulated gain of the signal in power writes as 



g(Q, z')dz' 



(4) 



with p = Pi(Q,0)/P s (Q,0). In the following, we set p = 1 for the sake of simplicity. In 
Eq. [4j g(Cl, z) = 2^P P sin[#(fi, z)\ is the local linear gain. Its calculation requires integrating 
Eq. 3c in order to evaluate 9(Q, z), which gives (under our assumptions) 

— B A VL 2 
6{Q, z) = [/3 2 ft 2 + 2 7 P p ]z + ^=[1 - cos(2%z/Z)} + 6({l, 0) (5) 

Z7T/Z 

Finally, by using a Fourier series expansion to calculate the sin[#(f2, z)\ term, we find that 
the local linear gain writes as 



' (6) 



g(Q, z) = 2 lPp J2 J, (|^) + 2^ - ^) 

Q=— oo \ / / L \ / 



Z + K q 



with K q = t^jz — g| + 9(Q,0). Thus, Eq. ^ indicates that the linear gain g(Q, z) at a 
fixed pulsation detuning Q can be interpreted as the sum of sine functions in z. These sine 
functions all have a zero average value except when their argument becomes independent 
of z. It only occurs at specific spectral components Vt (equal to the pulsation VLk in Eq. [2]) 
corresponding to solutions of the quasi-phase-matching relation [2j For these specific pul- 
sation detunings each term of the sum in Eq. Q leads to periodical amplification and 
deamplification phases along the DOF except for the uniform contribution corresponding to 
q = k. This last term therefore prevails over the other ones on the gain G(Q, z) for long 
enough propagation distances. Thus the linear gain of the k th spectral component mainly 
depends on this uniform term as long as the fiber exceeds a few modulation periods. It is 
then equal to 2jP p \ J k { 2 ^^ )|, by choosing K k = +| as initial condition. Note that this is 
analogous to the choice of maximizing the gain in MI in uniform fibers by setting the initial 
phase mismatch value to tt/2 
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FIG. 2: (a) Evolution of the gain of the first amplified frequency (k = 1) from Eq.[4j in red line with 
the contribution of Ji+Jo+J_i (average gain + oscillating terms), in blue solid line for J\ only, 
in blue dotted line for Jo only and in blue dashed line for J_i only with fijj™ = 2-ir x 2.63 x 10 12 
rad/s (Eq. [2J. The solid black line is calculated from the numerical integration of the original set 
of equations (Eqs. [3} with U s ™£ = 2ir x 2.93 x 10 12 rad/s. (b) Evolution of 6(z) from our analytical 
study (solid red line) and from numerical simulations with Eqs. [3] (solid black line). The green area 
corresponds to amplification and the red one to deamplification over one period. 

B. Physical interpretation 

To illustrate this process, we firstly focus on the first spectral component (k = 1). The 
solid black line in Fig. |2](a) shows the evolution of the maximum gain (at Q = f^ax ) obtained 
from numerical integration of the complete set of original Eqs. [3] Note that an excellent 
agreement is achieved with the numerical integration of the GNLSE (not shown here for the 
sake of clarity). The blue solid line in Fig. [2]^a) corresponds to the term of uniform gain 
(Bessel function Ji), the blue dotted and dashed lines correspond to the highest amplitude 
oscillating terms (Bessel functions Jo and J_i respectively in this case) and the red solid line 
corresponds to their sum. We limit our investigations to Jo and J_i because all other Bessel 
functions have much lower contributions in this example. A good agreement is obtained 
between the red solid curve from the analytical model, and the black one from numerical 
simulations, which confirms the validity of our assumptions and the accuracy of our method. 
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FIG. 3: (a) Evolution of the gain of the second amplified frequency (k = 2) from Eq. |4j in 
red line with the contribution of J-2+J-1+J0+J1+J2 (average gain + oscillating terms), in blue 
solid line for J 2 only, in blue dotted line for J\ only and in blue dashed line for J_i only with 
^max = 27r x 4.49 x 10 12 rad/s (Eq. [2]). The solid black line is calculated from the numerical 
integration of the original set of equations (Eqs. [3]) with fi^™" = 2-7rx4.61 x 10 12 rad/s. (b) Evolution 
of 9(z) from our analytical study (solid red line) and from numerical simulations with Eqs. [3] (solid 
black line). Green areas correspond to amplification and the red ones to deamplification over one 
period. 

In each modulation period, the amplification phase is characterized by < 9(Q, z) < tt and 
the deamplification one has — n < 6(VL,z) < as represented in Fig. [2]^b), the total phase 
shift being equal to 2ti per period. 

The dynamics of the second spectral component (k = 2) is shown Figure [3j The same 
reasoning as for the first one (k = 1) can be applied. The Ji term provides the average 
exponential gain (blue solid line) and additional oscillating terms provide the oscillating 
behavior of the overall gain. By taking the five terms with highest amplitude into account 
(from k = —2 to k = 2) a good agreement between Eq. [4] and numerical simulations from 
Eqs. [3] is achieved. Note that adding other higher-order terms does not significantly change 
analytical results (displayed in red curves). For the sake of clarity, only the two highest 
amplitude ones (J_i and J\) are represented in Fig. [3]^a) (in dashed and dotted lines respec- 
tively). In this case, there are two amplification and deamplification phases per period. The 
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evolution of the phase represented in Fig. |3[b) also shows a more complex evolution than 
for the first (k = 1) sideband and the total phase shift is now equal to 4n per period. The 
agreement between numerical simulations from Eqs. [3] (black lines) and our analytical result 
(red lines) is here excellent both for the evolution of the gain and of the phase. 

Figures [2] and [3] emphasize that the dispersion modulation enables to control the evolution 
of the relative phase of the waves so that the whole process can be seen as quasi-phase- 
matched, the variation of the relative phase over one period being equal to 2kir for the k th 
spectral component. Indeed, this relative phase 9(Q, z) would grow linearly in the absence of 



the modulation term = 0) in Eq. 3c [25] . This linear growth would lead to amplification 
and deamplification phases of same length and consequently the total accumulated gain 
would be negligible. 

From a more practical point of view, the frequency of the spectral component can 
be widely modified simply by changing the periodicity of the grating (as in a diffraction 
grating for the position of its different orders), while the gain (analogous to the diffraction 
efficiency in a specific order) can be modulated independently through the ratio /3^//3 2 - 
Note however that the deamplification phases of the signal along the fiber cannot be totally 
avoided but only reduced. This can be understood either by considering that they are 
due to the contribution of all oscillating terms of Eq. [6] which cannot be all suppressed 
simultaneously, or by considering that the dispersion grating enables the evolution of 9(Q, z) 
to deviate from a linear growth but not to limit its evolution in the [0; it] range required for 



a positive gain (see Eq. 3a). 



IV. EVOLUTION OF THE MI SPECTRUM WITH AVERAGE DISPERSION 

In order to further emphasize the accuracy of our simplified analytical model, we studied 
the evolution of the MI spectrum as a function of average dispersion value at the pump fre- 
quency both in normal and anomalous dispersion regimes. We chose the DOF parameters 
so that they match the ones of the fiber used hereafter. We took into account longitudinal 
variations of (3 2 according to Eq. [TJ and average values of 7, /3 3 and /3 4 instead of their longi- 
tudinal evolution. We checked numerically that this has negligible impact in our conditions. 
The DOF parameters are $ of 1(T 27 s 2 /m, ft = 6.8x KT 41 s 3 /m, & = 1.7 x 1(T 55 s 4 /m 
and 7 = 7 W _1 .km _1 at 1064 nm. The DOF is 120 m long and the modulation period is 
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FIG. 4: Top view of the gain spectrum (in logarithmic scale) obtained from (a),(c) GNLSE sim- 
ulations and from (b),(d) our analytical model (Eq. [4]) as a function of the average second-order 
dispersion. (a),(b) panels corresponds to average normal dispersion pumping and (c),(d) to anoma- 
lous average dispersion pumping. The zero frequency corresponds to the pump. 

Z = 10 m. Fiber attenuation and stimulated Raman scattering are neglected. 

We have firstly performed numerical simulations using Eqs. [4] by varying values of 
from -1.5 to 1.5xl0~ 27 s 2 /m. The maps displayed in Fig. |4](a) and (c) represent calculated 
numerical gain spectra (in logarithmic scale) for varying /3 2 values in normal and anomalous 
dispersion regions respectively. For the sake of clarity, only one half of the overall spectrum 
is displayed, the other half being perfectly symmetric with respect to the pump frequency. 
In parallel, we have plotted in Fig. |4](b) and (d) the same map using the analytical gain 
expression given by Eq. [4] with the same parameters as for the numerical simulations. There 
is a good qualitative agreement between both graphs, although quantitative agreement is 
not reached due to our approximations explained above. However, our model allows to 
reproduce the specific dynamical features observed in numerical simulations: 
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• Firstly, the detuning of each MI sideband from the pump decreases as /3 2 increases, as 
expected from Eq. [2} 

• Secondly, the gain strongly depends on /3 2 , which differs from classical MI spectrum 
observed in uniform fibers in the anomalous dispersion regime, in which the maximal 
parametric gain does not depend on /3 2 . Additionally, the maximal gain does not 
correspond to the same /3 2 value for each MI sideband. 

• Thirdly, we can identify specific /3 2 values for which one or several initially well-defined 
MI sideband are cancelled, i.e. their parametric gain vanishes. The frequency of 
cancelled sidebands increases with increasing /3 2 values, for each value of k. This 
a priori unexpected cancellation of parametric gain is in fact due to vanishing Jk 
functions in Eq [6] and will be detailed in the last section of this paper. 



V. TAILORING THE MI GAIN SPECTRUM 



The simple analytical approach presented in section |III| and further confirmed in section 
IV allows a better understanding of the complex dynamics of the process and it has allowed 
us to design experiments in which the multi-peak gain spectrum is tailored. To illustrate 
this, we focus our attention here to two striking examples. We chose either to completely 



cancel a given spectral component (as previously mentioned in section IV) or to maximize 
the gain of a sideband pair with regards to the others. To reach these goals, let us recall 
that, as detailed in the analytical model above, the linear gain of the k th spectral component 
can be approximated by 



g{Vt k ,z) = 2>yP p 



Jk. 



2-k/Z 



2 1 P P 



Jk 



IP2 



7T 



(7) 



Equation [7] indicates that the gain of the k th spectral component can be totally cancelled 
by simply finding the argument 77 = ^J/z = — ^ y ^ or wn i c h the Bessel function 

vanishes, 77 can be adjusted by controlling the modulation amplitude of dispersion, (3^, or 
the fiber period, Z, which both require manufacturing new DOF samples. But it can also 
be adjusted by controlling the average dispersion at the pump wavelength, fa, which can be 
done experimentally by simply tuning the pump wavelength. 
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A. Cancellation of spectral component 



To illustrate this, we fabricated a new DOF sample, labeled DOF^2 hereafter. DOF^2 
is 120 m-long and it has an average dispersion /3 2 of 10~ 27 s 2 /m, a modulation amplitude 
02 of 10 -27 s 2 /m at 1064 nm and its ZDW oscillates between 1064 nm and 1080 nm. The 
average third order dispersion term is 03 = 6.8 x 10 -41 s 3 /m, the average nonlinear coefficient 
is 7 = 7 W _1 .km _1 and the attenuation is a = 7.5 dB/km at 1064 nm. To illustrate the 
cancellation of parametric gain at specific frequencies, we choose for example to cancel the 
k = 6 sideband pair. In this case, we found that J 6 vanishes for a 02 of 5.8 x 10~ 28 s 2 /m 
with the parameters of DOF^2 given above and a pump power of 13 W. Circles in Fig. ga) 
shows the gain calculated with the above model (Eq. [7]) for each spectral component and 
for a fiber length of 120 m, while the solid line represents the output spectrum obtained 
from a numerical integration of the GNLSE with a pump power of 13 W. We average 50 
output spectra seeded by random initial conditions to account for the averaging performed 
during the experimental recording of a spectrum. These results firstly confirms the ability 
of our simplified model to correctly predict the maximal gain of each sideband and they 
also show that the k = 6 spectral component is indeed canceled. Experiments performed 
in DOF^2 by tuning the pump wavelength to 1067.5 nm (which is close to the required 02 
value of 5.8 x 10~ 28 s 2 /m) are displayed in Fig. gb). The overall shape of the experimental 
spectrum nicely matches the one obtained from theory and this measurement also confirms 
the cancellation of the 6 th peak. In all experiments presented in this section, the pump 
power was the only adjustable parameter. It had to be increased up to 24 W to observe 
the expected behaviors, which is higher than the power of 13 W used in simulations and 
in the model. This discrepancy in pump power is reasonable given the uncertainty on the 
evaluation of fiber properties (attenuation, dispersion and nonlinearity) and on the pump 
laser parameters (repetition rate, pulse duration, measurement of average power). 

B. Maximization of a single spectral component 

In order to further illustrate the possibility of tailoring the multi-peaks MI spectrum, 
we used Eq. [7] to find a configuration in which the k — 1 sideband is maximized, i.e. it 
experiences a much higher gain than any other ones. In this case, we simply need to find a 
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FIG. 5: Illustration of the cancellation of the k = 6 spectral component, (a) Maximal gain obtained 
from Eq. [7] (circles, right axis) and output spectrum simulated with the GNLSE (solid line, left 
axis), for fit = 5.8 x 1CF 28 s 2 /m and P p = 13 W. (b) Corresponding experiments performed in 
DOF#2 for a pump wavelength of 1067.5 nm and pump power of 24 W. 

02 value (and thus a value of the rj parameter), experimentally a pump wavelength, which 
maximizes the J\ Bessel function. Figure [6]^a) shows the gain calculated from Eq. [7] in 
circles, as well as the output spectrum obtained from numerical integration of the GNLSE 
for a @2 value of 3.87 x 1CT 28 s 2 /m. These results are again in excellent agreement, and 
they indeed show that the first sidelobe is favored since it has a 25 dB gain higher than 
all other ones. Experiments were performed by accordingly tuning the pump wavelength 
to 1071.5 nm. The output spectrum plotted in Fig. |6](b) shows that the power of the first 
sidelobe is 22 dB higher than other spectral components, in good agreement with theoretical 
predictions. 

It is also worth noting that in this case, the argument of the Ji Bessel function has been 
chosen so that the gain value of the first sideband calculated from our model corresponds 
to exactly the maximum value of 0.582 for J\. This means that the gain brought by the 
periodicity for this sideband approximately equals 0.582 x2 r jP p according to Eq. [7j It is 
less than a factor of 2 lower than the maximal gain expected from a classical MI process 
in the anomalous dispersion region in uniform fibers (which would approximately be equal 
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FIG. 6: Illustration of the maximization of the k = 1 spectral component, (a) Maximal gain 
obtained from Eq. [7] (circles, right axis) and simulated output spectrum (solid line, left axis), for 
@2 = 3.87 x 10 -28 s 2 /m and P p = 13 W. (b) Corresponding experiments performed in DOF^2 for 
a pump wavelength of 1071.5 nm and pump power of 24 W. 

to 2 r )P p ). This observation is all the more important given that no gain is expected in the 
normal dispersion region in uniform fibers (neglecting higher-order dispersion terms and 
higher-order fiber modes). 

Although the control of the overall MI spectrum shape requires a change of quasi-phase- 
matched frequencies, these examples demonstrate the possibility of harnessing the MI spec- 
trum thanks to the periodic dispersion landscape. A simultaneous control of both the 
spectral shape and sideband frequencies would still be possible by simultaneously adjusting 
$2 and 02 i which would however require manufacturing new fibers. 

VI. CONCLUSION 

Following our first experimental demonstration of MI in dispersion oscillating fibers p2] , 
we have investigated this process theoretically here. Starting from the well-known truncated 
three-wave model, we have derived an approximate analytical expression to predict the local 
parametric gain. This simplified model gives good agreement with numerical simulations 
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and experiments. It has also allowed to interpret the MI process in terms of a quasi-phase- 
matching mechanism due to the periodic nature of the fiber dispersion landscape. We 
have also used this model to emphasize the possibility of tailoring the MI spectrum, which 
has been confirmed experimentally by the cancellation or maximization of chosen spectral 
components. 

Dispersion oscillating photonic crystal fibers such as the ones reported here and in [T7j 
pave the way to a range of linear and nonlinear guided wave optical processes thanks to 
the longitudinal periodic modulation of their waveguiding properties. They should find 
applications in wavelength conversion, parametric amplification, generation of ultra-short 
pulse trains or soliton management among others. 
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